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STABILITY  or  A  PRISMATIC  BAR  UNDER  THE 
INFLUENCE  OF  A  COlfPRESSIONAL  STRAIN  WAVE 

In  this  paper  a  study  is  made  of  critical  strain  distributions  introduced 
into  a  bar  by  dynamical  loading  in  an  axial  direction.  Transverse  vibrations  of 
the  bar  which  may  result  from  this  type  of  loading  are  not  investigated.  It  is 
considered  that  they  are  really  Involved  in  an  extension  of  the  problem. 

The  theory  is  developed  in  some  detail  for  cases  of  bars  of  elastic  material 
having  Hookean  or  Non-Hookean  stress-strain  laws.  Also,  consideration  is  given 
to  the  case  of  the  bar  in  the  plastically  defonned  state. 

One  of  the  conclusions  from  the  study  is  that  the  proposed  method  for  deter¬ 
mining  critical  states  of  strain  suggests  a.  means  for  experimentally  investigating 
various  plastic  wave  theories.  Assuming  tentatively  a  given  theory  of  plastic 
strain  propagation  into  the  bar  the  critical  strain  distribution  may  be  calculated. 
Indications  of  plastic  buckling  may  then  be  investigated.  The  buckled  shapes  should 
also  provide  information  on  the  theory  of  wave  propagation. 

INTRODUCTION 

There  is  a  growing  interest  in  the  problem  of  structxures  subjected  to  dynamic 
loads  of  the  transient  type.  These  loads  are  sometimes  produced  by  collision  or 
explosion.  The  colunm  Is  an  Imnortant  structural  element  that  may  be  subjected  to 
this  kind  of  force.  The  case  of  periodic  axial  loads  of  the  steady-state  vibration 
type  has  been  considered  by  Stoker  and  Lubkin  (l) .  The  column  compressed  by  the 
moving  head  of  a  testing  machine  has  been  treated  by  Hoff  (2).  In  these  cases  the 
load  is  considered  to  vary  with  time  but  to  be  the  same  at  all  sections  of  the 
column  at  any  one  time.  In  other  words,  waves  of  coraoressional  strain  ere  not 
considered.  In  the  present  paper  a  study  will  be  made  of  the  compressional  strain 
waves. 
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DirFFRINTlAL  PTOATIONS  OF  INSTABILITY 


If  th«  stress  is  not  proportional  to  strain  we  may  define  the  nne-dimensional 
strain  in  the  manner  of  von  Karman  (3)t 


If  stress  is  proportional  to  strain  the  one-dimensional  wave  equation  may  be 
written  in  the  usual  manner: 


E. 


With  given  initial  end  boundary  conditions  for  a  bar  these  equations  may  be 
solved.  At  any  time  t  a  determinable  stress  distribution  exists  in  the 
bar. 

As  can  easily  be  shown  the  differential  equation  for  the  transverse  vibration 
of  a  longitudinally  loaded  bar  of  variable  stiffness  per  unit  length  may  be  written: 


Tt  ^  - 


o 


If  we  consider  \^lj  or  [^2j  and  [33  simultaneously  we  define  possible  trans¬ 
verse  displacements  induced  by  the  compressional  strain  wave.  In  order  to  obtain  a 
first  criterion  of  instability  it  is  proposed  to  neglect  the  acceleration  term  in 
^3^  and  determine  P  which  is  a  function  of  x  and  of  t  from  j^ljor  £2j.  The 
merit  of  these  assumptions  will  be  investigated  in  a  series  of  experiments  at  a 
later  date. 

The  D'Alembert  principle  can  be  used  and  the  P  consider  ad  as  reversed 

effective  forces  putting  the  bar  in  statical  eauilibrium  along  its  length.  This 
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principle  provides  a  means  of  determining  the  load  distribution  along  the  bar  or 
column  at  an7  time.  We  may  then  consider  a  force  at  the  wave  front,  a  force  at 


the  externally  loaded  end,  and  a  distribution  of  forces  between  these  two  points. 
The  force  corresponding  to  each  differential  length  is  ^AcJyc  •  This  theory 

will  be  illustrated  by  three  cases. 


CASE  I 

ELASTIC  CONDITION,  FOOKEAN  STRESS-STRAIN  LAW 
Assvune  a  straight  bar  initially  \jnstressed  end  with  certain  given  boundary  con¬ 
ditions.  The  material  of  the  bar  is  elastic  with  a  linear  stress-strain  law.  A 
compressive  load  as  a  known  function  of  time  is  suddenly  applied  at  one  end.  If 
numerical  methods  of  solution  are  to  be  avoided,  and  they  would  be  quite  clumsy  in 
a  problem  of  this  t3rpe,  analytical  solutions  must  be  obtained.  Many  of  these  can 
be  found  in  the  literature.  The  simplest  will  now  be  considered  (4)» 

In  Fig.  1  a  simply  supported  bar  of  length  1  is  shown  with  a  force  P  applied 
at  one  end  and  the  front  of  the  strain  wave  progressed  to  a  distance  Xq  from  that 
end.  Assvune  that  the  force  P  is  caused  by  a  large  mass  M  moving  with  constant 
velocity  Vq.  Then  the  strain  will  be  constant  and  of  magnitude  given  by 


where  c  is  the  constant  velocity  with  which  the  strain  moves  into  the  bar  (5). 

The  principle  of  D'Alembert  will  now  enable  us  to  translate  this  problem  into 
an  equivalent  static  one  for  any  given  time  t.  The  corresponding  loading  for  this 
case  is  then  shown  in  Fig.  2. 
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P  =  AEK 


In  region  ©  t^ere  If  no  Ftrein  since  the  strein  T/sve  has  not  progressed  into 

that  oortion  of  the  bar  6?  yet.  In  Fegion  there  Is  constant  strain  of  nagnitude 

at  eac’"  noint  and  load?  of  magnitude each  of  its  extremities. 

^  C 

The  two  differential  eo^oations  for  bending  of  the  bar  are: 


Region  (T) 


El  p  =  a 

ax'*-  ^ 


Pegiondp . 


Solution  of  [*+]  is: 

w,  =--  A  ,  ^  X  -b  X) 

and  solution  of  [3Q  is 

w\/  —  J\  -f-  S  ~i~  C.  5  cw  A  X  “4-  A  X 

'  z.  ^  i~  jL 

The  Boundary  and  continuity  conditions  are: 


C^l 
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The  vanishing  of  the  detemiinant  of  the  system  of  the  8  linear  eqiiatlons  in 
the  8  constants  A^,  Bi,  Ci,  Di,  A2,  C2,  I?2  the  instability  equationi 


Co-^  A  ^2  ■+■  S  ^  '^2  ~  ^ 


rsi 


For  convenience  we  may  put 


then  equation  [s^  becomes: 

-  f  )  p  I®  ?  +  ^  f  -  o 


I’J 


67  means  of  aquation  F^J  the  following  table  of  critical  loads  is  determined. 


f 

TABLE  I. 

A 

_  P 

A  -  ^ 

0. 

00 

C^nO 

.1 

16.35 

268.3 

.2 

8.55 

73.0 

.3 

6.01 

36.1 

.4 

4.76 

22.7 

.5 

4.06 

16.4 

.6 

3.62 

13.1 

.7 

3.36 

11.2 

.8 

3.21 

10.3 

.9 

3.15 

9.9 

1.0 

3.U 

9.9 

1?^ 


[10] 


It  is  readily  seen  that  a  relationship  has  been  established  between  the  fractional 
length  factor  f  and  the  critical  load  given  by  [loj .  If  reflections  of  the 

strain  wave  occtir  before  the  critical  stress  is  attained,  the  process  must  be 
repeated  with  appropriate  strain  distributions  and  another  table  of  f  values. 


CASE  II 

ELASTIC  CONDITION,  NOM-HDOKEAN  STRESS- 
STRAIN  LAW 

In  thecase  of  bars  made  of  non-Hookean  materials  the  stress-strain  law  is  non- 
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I 


Fig.  3.  Stress-strain  law  as  two  straight  line  segments 


Fig,  4.  Strain  distribution  in  bar 


linear,  but  the  materials  may  still  he  elastic.  A  simple  example  of  such  a  stress- 
strain  relation  is  shown  in  Fig.  3* 

As  an  example  of  this  case  suppose  a  force  is  applied  at  one  end  of  a  straight 
bar  so  that  the  load  is  constant  in  time  and  produces  a  stress,  at  the  end  of  the 
bar,  which  lies  between  and  in  Fig.  3*  In  this  case  a  wave  of  stress 

of  magnitude  oropagates  at  a  known  speed  into  the  bar.  This  wave  is  followed 

by  a  slower  one  of  magnitude  as  shown  in  Fig.  4«  In  Fig.  5  the  equivalent 

static  loading  at  some  time  t  is  shown. 

In  region  the  load  is  a  force  (Pl^“P2)  in  region  the  load  is  and 
in  region  (T)  there  is  no  load  as  the  wave  front  has  not  penetrated  this  region  as 
yet.  This  strain  propagation  satisfies  the  differential  equation  [[l]  .  The  dif¬ 
ferential  equations  of  lateral  bending  are; 


E  I 


Region 


El 


o 


Region 


El 


=  .O 


Region 
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The  solutions  are: 


=  4 ,  ■(-  S,  X  c  sov,  >  .X  +  ^ 

-*■  sv^  -I- 

Wj  =  x\  x^+  X  +  i:::^ 


■■■  C 


The  boundery  and  continuity  conditions  are: 


and 


■) 


-  -  -  U<^1 


O'J 


The  venlshing  of  the  determinant  of  the  eystem  of  12  linear  homogeneous  equations 
in  the  12  constants  gives  the  instability  equation: 

—  O  -  -  -  -  - - II'2^ 


^e  may  divide  equation 

ratios  'cr-' 

^  1 

and 

CJ^  / 

-  ^  - 

2 

^  nc 

and  determine  the  roots  for  various 
As  an  example  suppose: 


Then  we  may  write  as  follows: 


C  &  - 

c  t  =  ^, 

a-  f 

II 

=/^ 
/  U  ^ 

-  8  - 


^ -/ 

/  3-  J  J  '  i 


Then 


&T-J< 

-  ^-  /  -  Z7S 


Then  let  1 


where 


/  ^ 


is  R  fractional  parameter. 


4- 


For  the  assumed  ratios 


Cj/  / 
and  '  . 


we  have 


P,  +  P  ^ 

A  "  ^ 

'-  f 


^>//-/r/// 


9 


Fig.  5*  Equivalent  static  loading 


€ 


Fig.  6.  Non-linear  stress-strain  curve  of  strain 
bardenable  material 


Obviously  the  nexlnum  slloTOble  In  this  formulation  is  given  by 

Then  re  ney  construct  c  table  of  values  relating  critical  loads  to  lengths. 


TABLE  II 


.1 

7.48 

55.90 

14.96 

223.6 

.2 

3.94 

15.40 

7.88 

61.6 

.3 

2.78 

7.70 

5.56 

30.8 

.4 

2.24 

5.00 

4.48 

20.0 

.5 

1.94 

3.75 

3.88 

15.0 

.6 

1.76 

3.09 

3.52 

12.4 

.7 

1.65 

2.71 

3.30 

10.8 

.7071 

1.64 

2.68 

3.28 

10.7 

Obviously  for  any  given  force  P, 

the  critical  : 

fractional  length 

^  is 

determined 

.  The 

shape  of  the  bent  bar 

is  given  by  equatipns  j_v  J 

.  Again 

in  this  case  if  reflections  of  the  strain  wave  occur  before  the  critical  stress 


is  attained  the  process  must  be  repeated  with  appropriate  strain  distributions 
and  smother  table  of  ^  values. 

CASE  III 

PLASTIC  CONDITION 

The  differential  equation  of  the  stjrain  in  bar  in  the  plastic  region 
exhibiting  unloading  phenomenon  may  be  written  as  noted  by  T.  von  Kannan 


^  _  O  .PV 


where  ^  ^  is  the  slope  of  the  stress  strain  curve  and  is  the  density « 
It  is  claimed  that  the  velocity  of  such  a  plastic  wave  is 


-  /-L  <^0" 

y  r 


and  the  static  engineering  stress-strain  curve  for  the  material  may  be  used 
to  determine  ,?2r  for  any  given  ^  . 

As  an  example  suppose  that  a  bar  shown  in  Pig*  1  is  subjected  to  a  con¬ 
stant  stress  St  one  end.  Then  it  is  asserted  that  a  plastic  wave  of 

intensity  CT^  progresses  into  the  bar  at  velocity 


rr^ 


In  order  to  determine  instability  in  this  case  for  a  material  exhibiting 
strain-hardening,  we  may  use  the  method  of  Elngesser  (6)  to  take  account  of 
unloading. 

Recalling  the  method  of  Engesser  an  effective  modulus  is  determined  for 
the  bar  section  as  in  Fig.  7. 

.  K  =  E,K,  .  >ij 

ySTS  Z.  j^*  ^ 


t 

Z>  ^ 


30 
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which  is  the  force 
equation 


i 


K.  -  K 


.E,  2J1 


a.  -  (M 


M  is  the  moment  at  section. 

If  b  Is  the  width  of  the  section 

„  =  Ail .  4 

12  f 

piutting  pT"  _  4- 


arjid  M  -  Ej.  *  ^ 
yO** 


If  we  now  replace  E  by  Ej.  we  may  revert  to  Case  I  for  determining  the 
l.ength  of  loading  section  for  any  given  applied  end  stress  CT^  which  is 
co'i'tlcal. 

DISCUSSION 

The  theory  presented  herein  is  proposed  as  a  first  criterion  of  instability 
o«f  a  prismatic  bar  with  a  transient  force  applied  at  its  end.  Although  a  complete 
md  accurate  equation  describing  strain  wave  propagation  in  a  bar  must  include 
Lateral  contraction  as  well  as  other  effects  in  the  bar  it  is  considered  that  the 
ssimplest  one-dimensional  wave  theory  serves  the  purpose  of  the  present  exposition. 
For  a  similar  reason,  the  transverse  inertia  forces  related  to  bending  are  neglected 
Un  order  to  get  an  approximate  criterion  of  instability.  It  is  considered  that 
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The  merit  of  these  asstimptlons  should  be  experimentally  Investigated,  because  very 
few  data  on  the  compression  Impact  loading  of  bars  occur  In  the  literature.  Some 
experiments  by  D.  S.  Clark  (7)  are  quite  sviggestlve  of  what  may  be  expected  In  the 
way  of  sidewise  plastic  deformation  of  relatively  long  bars  of  lead  subjected  to 
dynamic  loads  on  one  end.  A  series  of  experiments  will  be  conducted  now  by  the 
author  with  the  specific  purpose  of  investigating  sidewise  deformation  of  bars 
subjected  to  Impact  loading  on  one  end. 
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